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Abstract 

In this paper, we analyse the joint impact of pathloss, shadowing and fast fading on wireless networks. Taking into account 
the pathloss and the shadowing, we first express the SINR distribution of a mobile located at a given distance from its serving 
base-station (BS). The moments of this distribution are easily computed, using the Fenton- Wilkinson method, and a fluid model 
that considers the cellular network as a continuum of BS. Then considering the joint impact of pathloss, shadowing and fast 
fading, we derive an easily computable outage probability formula, for a mobile located at any distance from its serving BS. 
We validate our approach by comparing all results to Monte Carlo simulations performed in a traditional hexagonal network. 
Indeed, we establish that the results given by the formula are close to the ones given by Monte Carlo simulations. The proposed 
framework is a powerful tool to study performances of cellular networks e.g. OFDMA systems (WiMAX, LTE). 



I. Introduction 

The characterization of interference allows to analyze cellular networks performance. In cellular networks, an important 
parameter for this characterization is the SINR (signal to interference plus noise ratio) at any point of the cell. Indeed, it allows 
the derivation of performance evaluation (outage probabilities, capacity,...). We focus on downlink throughout this paper, since 
that direction is often the limited link in term of capacity. Nevertheless, our approach can easily be extended to the uplink. 

The issue of expressing outage probability in cellular networks has been extensively addressed in the literature. For this 
study, there are two possible assumptions: (1) considering only the shadowing effect, (2) considering both shadowing and fast 
fading effects. In the former case, authors mainly face the problem of expressing the distribution of the sum of log-normally 
random variables; several classical methods can be applied to solve this issue (see e.g. [14] [3]). In the latter case, formulas 
usually consist in many infinite integrals, which are uneasy to handle in a practical way (see e.g. [15]). In both cases, outage 
probability is always an explicit function of all distances from the user to interferers. 

As the need for easy-to-use formulas for outage probability is clear, approximations need to be done. Working on the 
uplink, [13] derived the distribution function of a ratio of path-losses with shadowing, which is essential for the evaluation of 
external interference. For that, authors approximate the hexagonal cell with a disk of same area. Authors of [16] assume perfect 
power control on the uplink, while neglecting fast fading. On the downlink, Chan and Hanly [11] precisely approximate the 
distribution of the other-cell interference. They however provide formulas that are difficult to handle in practice and do not 
consider fast fading. Immovilli and Merani [17] take into account both channel effects and make several assumptions in order 
to obtain simplified formulas. In particular, they approximate interference by its mean value. Outage probability is however 
an explicit function of all distances from receiver to every interferes Zorzi [8], proposes a formula available but for packet 
radio networks rather than for cellular systems. Authors of [9] provide some interesting characterizations and upper bounds of 
the outage probability but neglects the slowly varying pathgains. In [10], authors consider both shadowing and fast fading but 
assume a single interferes 

In this paper, to calculate the SINR distribution of a mobile located at any distance from its serving base station (BS), 
we consider the useful power he received and the interferences coming from all the BS stations of the network. We jointly 
take into account pathloss, shadowing and fast fading and derive a simple and easily computable outage probability formula. 
We validate that formula by comparisons of the figures obtained with it to the ones obtained by Monte Carlo simulations. 
At last, we rely on a recently proposed fluid model [5] [6] in order to express the outage probability as a simple analytical 
expression depending on the distance to the serving BS. Such an expression allows further integrations much more easily than 
with existing formulas. The formula can be used to analyze the coverage of a BS. 

In the next section, we derive outage probabilities, while considering first only pathloss and shadowing and then pathloss, 



shadowing and fast fading jointly. The computation is based on the fluid model (section IIIi. In section IV we validate our 



approach and compare analytical expressions with results obtained through Monte Carlo simulations. 



II. Outage Probabilities 
We consider N interfering base station (BS), a mobile u and its serving BSq. 
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A. Pathloss and Shadowing impact 

1) Propagation: Considering the power Pj transmitted by the BS j, the power pj iU received by a mobile u can be written: 

Pi,u = I'.hrJ),.,,. (1) 

where Yj^ u = 10 i° represents the shadowing effect. The term Yj tU is a lognormal random variable characterizing the random 
variations of the received power around a mean value. £j tU is a Normal distributed random variable (RV), with zero mean 
and standard deviation, a, comprised between and 10 dB. The term PjKrJ^ where K is a constant, represents the mean 
value of the received power at distance r$ u from the transmitter (BSj). The probability density function (PDF) of this slowly 
varying received power is given by 

, . 1 /ln(s) — am\ 
PY(s) = r exp - -^L (2) 

where a = m = ^ \n(KPjr~^) is the (logarithmic) received mean power expressed in decibels (dB), which is related to 
the path loss and a is the (logarithmic) standard deviation of the mean received signal due to the shadowing. 

2) SINR calculation: Considering the useful power Pq transmitted by the BSq, the useful power pq :U received by a mobile 
u belonging to BSq can be written: 

Po,u = P Kr^ u Y Q<u . (3) 

For the sake of simplicity, we now drop index u and set r^ :U = r. The interferences received by u coming from all the other 
base stations of the network are expressed by: 

N 

y^i-,K ■ ,■/>)). (4) 

The SINR at user u is given by: 

If now, thermal noise is negligible and BS have identical transmitting powers, the SINR can be written 7 = 1/Yf with 

N 

3) Distribution ofYf. The factor Yf is defined for any mobile u and it is location dependent. The numerator of this factor 
is a sum of log-normally distributed RV, which can be approximated by a log-normally distributed RV [3]. The denominator 
of the factor is a log-normally distributed RV. Yf can thus be approximated by a log-normal RV. 

Using the Fenton- Wilkinson [1] method, we can calculate the mean and standard deviation, rrif and Sf, of factor Yf, for 
any mobile at the distance r from its serving BS, BSq. We first calculate the mean and the variance of a sum of log-normal 
RV (appendix 1). We afterwards apply the result to the sum of N log-normal identically distributed RV (appendix 2): 

m f = -]n(yf(r,r))H(r,a)), (7) 

s) = 2(a 2 - 4 In #0,(7)), (8) 



1 = „ „ -r,„ . 71 ■ ( 5 ) 



where 



H{r,o) = e a2 ° 2 1 2 (G{r, V ){e a2 ° 2 + % (9) 



-2v 



G(r,rj)= , J J 2 , (10) 



yf(r,v) = ^r- (id 



From Eq. |6]and 11 we notice that yf(r,rj) represents the factor Yf without shadowing. 
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4) Outage probability: The outage probability is denned as the probability for the SINR 7 to be lower than a threshold 
value 8 and can be expressed as: 

P( 7 <<$) =l-¥(p >S( Pext + N th )) (12) 



So we have: 



»( 7 <«5) = l-¥^>Y f (m f ,8 f ) 

= i-p(ioio glo (i)>ioio glo (y / )) 



Q 



(13) 
(14) 

(15) 



where Q is the error function: Q(u) = \ er f c {^) 

The outage probability for a mobile located at a distance r from its serving BS, taking into account shadowing can be 
written: 

"iOlogioU) -m f ~ 



D ( 7 < S) = Q 



(16) 



B. Pathloss Shadowing and Fast fading impact 

1 ) Propagation: The power received by a mobile u depends on the radio channel state and varies with time due to fading 
effects (shadowing and multi-path reflections). Let Pj be the transmission power of a base station j and Xj iU the fast fading 
for mobile u, the power p^ u received by a mobile u can be written: 



Pj,u 



P Kr~ v X- Y 

3 j,u y ^3-M L J,ui 



where X^ u is a RV representing the Rayleigh fading effects, whose pdf is p x (x) — e x . 
2) SINR calculation: The interference power received by a mobile u can be written: 



Pext,u — PjK r j,uXj,vYj,u- 



(17) 



(18) 



Considering that the thermal noise is negligible and that all BS transmit with the same power P , we can write (dropping 
index u): 

7= (19) 
7 22Ur/>X jYj (19) 



3) Outage probability: As previously, the outage probability is expressed as: 



N 



»( 7 < S) = 1 - ¥(r^X Q Y > S(J2 rpXjYj)). 



(20) 



The interference power received by a mobile u due to fast fading effects varies with time. As a consequence, the fast 
fading can increase or decrease the power received by u. We consider that the increase of interfering power due to fast fading 
coming from some base stations are compensated by decrease of interfering powers coming from other base stations. As a 
consequence, for the interfering power, only the slow fading effect has an impact on the SINR. So, we assume that Xj can be 
neglected V j 7^ (this assumption will be validated by simulations in the next section): F(r~ v XoY > <5(X)^=i r J n XjYj)) w 
F{r-^X Y >S(Ef =1 rpY j )). 



So we have: 



1 N 

»( 7 <«5) = l-F(X Y >5—(£rpY i )) 



3=1 



l-¥(X (t)>SY f ) 



= 1 
= 1 



¥(x > 8Y f )px{x)dx) 
P(l01og 10 (-) > lOlog^r^e-^. 
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As a consequence, the outage probability for a mobile located at a distance r from its serving BS, taking into account both 
shadowing and fast fading can be written: 



P(7 < S) = 

where Q is the error function: Q(u) = \ er f c {^) 



IQlogio(f) -m f 
s f 



c dx, 



(21) 



Remark 



We notice that Eq. 21 characterizes the cumulative density function (CDF) of 7. That formula also allows to analyze the 
coverage of a BS for a given outage probability. 



III. Analytical fluid model 

To calculate the outage probability we need to express analytically the factors mt and St. In this aim we use the fluid 
model approach developed in [7]. The fluid model and the traditional hexagonal model are two simplifications of the reality. 
As shown in [7], none is a priori better than the other but the latter is widely used, especially for dimensioning purposes. The 
key modelling step of the model consists in replacing a given fixed finite number of transmitters (base stations or mobiles) by 
an equivalent continuum of transmitters which are distributed according to some distribution function. We consider a traffic 
characterized by a mobile density and a network by a base station density pss- For a homogeneous network, the downlink 
only depends on the distance r between the BS and the mobile. We denote it yf(r). From [7], we have 



parameter yt (Eq. 



11 



Vf{r) 



V-2 



[(2R C - r) 2 -" - (R nw - r) 



2-7,1 



(22) 



where 2R C represents the distance between two base stations and R nw the size of the network. We notice that the factor G 



(Eq. 10 1 can be rewritten (generalizing Eq. 22 dropping r, and expressing y as a function of 77) as: 

Vf&V) 



Vf(v) 



(23) 



The higher the standard deviation of shadowing, the higher the mean and standard deviation (see hereafter Eq. 24 and 25 1 



This increase is however bounded as we can see hereafter. In a realistic network, a is generally comprised between 4 and 8 
dB. 

The Eq. [7] means that the effect of the environment of any mobile of a cell, on the yj factor, is characterized by a function 
H(r,a). This function takes into account the shadowing and a topological G factor which depends on the position of the 
mobile and the characteristics of the network as 

• the exponential pathloss parameter 77, which can vary with the topography and more generally with the geographical 
environment as urban or country, micro or macro cells. 

• the base stations positions and number. 

From Eq. 10 we have < G(r, 77) < 1 whatever r and r\. From Eq. [7] we can express: 



Vf(r,v) <m f < 



Vf( r , v) 
G{r,ri)* 



and 



s) < 2a 2 . 



(24) 



(25) 



For low standard deviations (less then 4 or 5 dB), we can easily deduce from Eq. [7] a low dependency of the mean value 

of y with a, and the total standard deviation Sf is very close to a. These expressions show another interesting result. We 

22 

established that < G < 1. Without the topological factor G (G = 0), we have H = e a a so that m/ increases indefinitly. If 
G = 1, then H = 1 and rrif = - ln(y^). For each given value of r\, characterizing a given type of cells or environment, the 
topological factor G thus compensates the shadowing effects. 



IV. Performance evaluation 



In this section, we compare the figures obtained with analytical expressions (Eq. 21 1 to those obtained by Monte Carlo 
simulations. We also validate our approach and show the limit of our formulas. 



Fig. 1. Hexagonal network and main parameters of the study. 




Fig. 2. Variations of mj [dB] with a [dB] for r\ = 3. 



A. Monte Carlo simulator 

The simulator assumes an homogeneous hexagonal network made of several rings around a central cell. Figure [T] shows an 
example of such a network with the main parameters involved in the study. 

The simulation consists in computing the SINR for any point u of the central cell. This computation can be done independently 
of the number of MS in the cell and of the BS output power (because noise is supposed to be negligible both in simulations 
and analytical study). At each snapshot, shadowing and fast fading RV are independently drawn between the MS and the 
serving BS and between MS and interfering BS. SINR samples at a given distance from the central BS are recorded in order 
to compute the outage probability. 



The validation consists in comparing results given by simulations to the ones given by the Eq. 21 



B. Results 

Figure [2] shows that in presence of shadowing, the value of mj increases. However, we observe that this increase reaches 
a maximum value for a = 12dB. It is due to the effect of the topology factor G which compensates the effect of shadowing. 
To validate our analytical approach, we present results for different pathloss factors r\ =3 and 4, different standard deviations 
a = 3 and 6 dB, and for mobiles located at the edge of the cell r — R c and the center of the cell r — i? c /2. 

Figure [3] and [4] give two examples of the kind of results we are able to obtain instantaneously thanks to the formulas derived 
in this paper for the path-loss exponents 77 =3 and 4. These curves represent the outage probability for a mobile located at 
the edge of the cell. Figure [5] gives outage probability considering a path-loss exponent 7/ =3, for a mobile at a distance of 
r = R c /2. Figures [3] [4] and [5Jshow that our analytical approach gives results (solid blue curves) very close to the ones obtained 
with simulations (dotted red curves). 
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SINR (dB) 



Fig. 3. Outage probability for a mobile located at the edge of the cell: r = R c ; comparison of the fluid model (solid blue curves) with simulations (dotted 
red curves) on an hexagonal network with r) = 3. 
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Fig. 4. Outage probability for a mobile located at the edge of the cell: r = R c ; comparison of the fluid model (solid blue curves) with simulations (dotted 
red curves) on an hexagonal network with r) = 4. 




SINR (dB) 



Fig. 5. Outage probability for a mobile located at r = 0.5i? c ; comparison of the fluid model (solid blue curves) with simulations (dotted red curves) on an 
hexagonal network with r\ = 3 and <r = 3dB. 



Figures [3] |4][5] and [6] show the specific impact of shadowing alone (Figure [6]l and both shadowing and fast fading (Figures 
[3] [4] |5]l on the outage probability. In particular the curve a = dB (Figure [6} shows that, without shadowing and fast fading, 
the SINR received at the edge of the cell reaches -4dB. Since there is no random variations of the propagation, this value is 
deterministic: the outage is 0% until a SINR of -4dB and 100% for higher values of SINR. When the shadowing is taken into 
account (but not the fast fading), we observe for a = 2>dB, a SINR of -8 dB for an outage probability of 10%. And when 
the shadowing and the fast fading are both taken into account (Figure pj, we observe for a — 3dB, a SINR of -16 dB for an 
outage probability of 10%. 

With these results, an operator would be able to analyze the outage, and to admit or reject new connections according to each 



-40 -30 -20 -10 10 20 30 40 

SINR (dB) 



Fig. 6. Outage probability for a mobile located at r = R c ; shadowing impact (without fast fading) for r) = 3. 

specific environment of the entering MS. Indeed, the SINR depends on the mobile location and on the environment (Eq. H), 
i.e. the shadowing, the exponential pathloss parameter r\ (which depends on different characteristics and particularly the cell 
dimensions and the type of environment e.g. urban or country), the number and positions of the base stations. 

C. Applications 

Hereafter we present some other examples of possible applications of the analytical approach we proposed. 

1) Coverage: The coverage of a BS can be defined as the area where the SINR received by a mobile reaches a given value 
with a given outage probability. Figure 4 shows, for r/= 4, that at the edge of the cell (r= i?, c ), the SINR reaches the value of 
-15 dB with a probability of 10%, when a = 3dB. With the same outage probability (10%), the SINR only reaches -20 dB 
when a = 6dB. 

2) Capacity: Since the capacity and the throughput D can be expressed as a function of the SINR, for example using the 
Shannon expression D cx log 2 (l + SINR), the analytical method we developed can be used to analyze this performance 
indicator instantaneously. 

V. Conclusion 

In this paper, we establish a simple formula of the outage probability in cellular networks, while considering pathloss, 
shadowing and fast fading. Taking into account pathloss and shadowing, we first express the inverse of the SINR of a mobile 
located at a given distance of its serving BS as a lognormal random variable. We give close form formulas of its mean m/ and 
standard deviation Sf based on the analytical fluid model. It allows us to analytically express the outage probability at a given 
distance of the serving BS. We then consider both pathloss, shadowing and fast fading and give an analytical expression of the 
outage probability at a given distance of the serving BS. The fluid model allows us to obtain formulas that only depend on the 
distance to the serving BS. The analytical model that we propose is validated by comparisons with Monte Carlo simulations. 
The formulas derived in this paper allow to obtain performances results instantaneously. 

The proposed framework is a powerful tool to study performances of cellular networks and to design fine algorithms taking 
into account the distance to the serving BS, shadowing and fast fading. It can particularly easily be used to study frequency 
reuse schemes in OFDMA systems. 

Appendix 1 

Each BS transmits a power Pj = P so the power received by a mobile is characterized by a lognormal distribution Xj as 
ln(Xj) tx N(amj,a 2 a-j) and we can write rrij — ^^(PjrJ 71 ) (to simplify the calculation we consider K=l). So the total 
power received by a mobile is a lognormal RV X characterized by its mean and variance ln(X) cx N(am, a 2 af) and we can 

write: am = In ^X^^Lo j=tb e^ lnPj ^ r i lnr i + ~r JL ^\ — 2-^*- so we have since Pj = P whatever the base station j, considering 

identical Oj = o: 

22 ^ 22 

am = (lnP+^)+ln( £ g^ 1 ^)-^ 1 (26) 

j=0,j^b 

We can express the mean interference power I ext received by a mobile as: 

22 •'v 22 

ln(l^ t ) = (\nP+ a ^-)+ln( £ r j") - ^ (27) 
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The variance a 2 a 2 of the sum of interferences is written as 

a 2 CT 2 =ln | ^_ ;| | (28) 




Introducing 

^hi^ll (29 ) 

{^3=0,j^b'j J 

the mean value of the total interference received by a mobile is given by 

/. _2_2 . . . n-i/2 



and 



I ex t = P E ((e^ 2 - l)G(r7) + l) (30) 

a 2 a 2 t = In ((e^ - 1)0(7?) + l) + In (e^ 2 ) (31) 



and finally, denoting: 

we have 

In dB, we can express 



Appendix 2 

We denote I int the power received by a mobile coming from its serving BS. Since the ratio of two lognormal RV's is also 
a lognormal RV, we can write the following mean and logarithmic variance: 

m f = = (32) 

lint 

and thus, considering that: Pb = P, we can write, dropping the index b: 

■ t ■ ^ 22/ , , \ — 1/2 

m f = (V CT - l)G(v) + l) (33) 

H(a) = (V V2 - l)G(r?) + l) ^ (34) 

m f = 2//( r7 ) i ^( cr ) (35) 

^f d B= l -HVf{ri)H(<r)) (36) 
In a analogue analysis, the standard deviation is given by a 2 s"j — a 2 a 2 + a 2 a 2 so we have 

a 2 s) = 2{a 2 (j 2 - ln(H(a))) (37) 

References 

[I] L.Fenton "The sum of lognormal probability distributions in scatter transmission system", IEEE (IRE) Transactions on Communications, CS-8, 1960 
[2] S. Shwartz and Y.S. Yeh, On the distributions functions and moments of power sums with lognormal components, Bell Syst, Tech J, vol 61, pp 1441 

1462, Sept 1982 

[3] G.L. Stuber, principles of Mobiles Communications, Norwell , MA. Kluwer, 1996 

[4] Adnan A. Abu-Dayya and Norman C. Beaulieu, Outage Probabilities in the Presence of Correlated Lognormal Interferes, IEEE Transactions on Vehicular 

Technology, Vol. 43, N 1, February 1994 
[5] J.-M. Kelif, Admission Control on Fluid CDMA Networks, Proc. of WiOpt, Apr. 2006. 

[6] J-M. Kelif and E. Altman, Downlink Fluid Model of CDMA Networks, Proc. of IEEE VTC Spring, May 2005. 

[7] J-M Kelif, M. Coupechoux and P. Godlewski, Spatial Outage Probability for Cellular Networks, Proc. of IEEE Globecom Washington, Nov. 2007 
[8] M. Zorzi and S. Pupolin, Outage Probability in Multiple Access Packet Radio Networks in the Presence of Fading, IEEE Trans, on Vehicular Technology, 
Vol. 43, No. 3, Aug. 1994. 

[9] J. Papandriopoulos, J. Evans, and S. Dey, Outage-Based Optimal Power Control for Generalized Multiuser Fading Channels, IEEE Trans, on 

Communications, Vol. 54, No. 4, Apr. 2006. 
[10] A. Williamson and J. Parsons, Outage Probability in a Mobile Radio System Subject To Fading and Shadowing, IEEE Electronics Letters, Vol. 21, No. 

14, Jul. 1985. 

[II] C. C. Chan and Hanly, Calculating the Outage Probabability in CDMA Network with Spatial Poisson Traffic, IEEE Trans, on Vehicular Technology, 
Vol. 50, No. 1, Jan. 2001. 

[12] S. E. Elayoubi and T. Chahed, Admission Control in the Downlink of WCDMA/UMTS, Lect. notes comput. sci., Springer, 2005. 
[13] J. S. Evans and D. Everitt, Effective Bandwidth-Based Admission Control for Multiservice CDMA Cellular Networks, IEEE Trans, on Vehicular 
Technology, Vol. 48, No. 1, Jan. 1999. 

[14] M. Pratesi, F. Santucci, M. Ruggieri, Outage Analysis in Mobile Radio Systems with Generically Correlated Log-Normal Interferers, IEEE Trans, on 
Communications, Vol. 48, No. 3, Mar. 2000. 



9 



[15] J.-P. M. G. Linnartz, Exact Analysis of the Outage Probability in Multiple-User Mobile Radio, IEEE Trans, on Communications, Vol. 40, No. 1, Jan. 
1992. 

[16] V. V. Veeravalli and A. Sendonaris, The Coverage-Capacity Tradeoff in Cellular CDMA Systems, IEEE Trans, on Vehicular Technology, Vol. 48, No. 
5, Sept. 1999. 

[17] G. Immovilli and M. L. Merani, Simplified Evaluation of Outage Probability for Cellular Mobile Radio Systems, IEEE Electronic Letters, Vol. 27, No. 
15, Jul. 1991. 



